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I. INTRODUCTION 

A point vortex system is a model of continuous two- 
dimensional (2D) inviscid fluid dynamics. In fact, a point 
vortex is an idealized vortex of real vortex in 2D fluids. 
However, the point vortex system carries some important 
properties of continuous 2D inviscid fluid dynamics]!! 2 1 . 
It is formulated as a Hamilton dynamical system [3j, |4|. 
Therefore, a standard statistical mechanical theory was 
developed by OnsagerQ. In his formulation, states of a 
point vortex system are classified into two categories: one 
is positive temperature states, and the other is negative 
temperature states. In the positive temperature states, 
point vortices are distributed in a scattered way. How- 
ever, in the negative temperature states, huge vortices 
are formed in time-evolution. The negative temperature 
states are phenomena observed in earth's surface. That 
is a formation of huge vortices, such as hurricanes and ty- 
phoons. Onsager's ideas are recently reviewed by Eyink 
and SreenivasanQ. Since Onsager, researchers consid- 
ered the equilibrium statistical mechanics of point vor- 
tex systems. The main purpose of their studies is to 
construct the equilibrium states of point vortex system 
0, B i, EE EL El- The Poisson-Boltzmann equation is 
used to obtain the equilibrium states. For some cases, 
the equilibrium states are analytically obtained [T^]. 

The next stage of the statistical mechanics of point vor- 
tex systems focuses on nonequilibrium properties. Stan- 
dard nonequilibrium statistical mechanics goes toward 
kinetic theory. Two authors have studied kinetic the- 



ory of point vortex system, i.e., by Marmanis 14,] and by 
Chavanis (To|. Marmanis considered a gas of binary pairs 
of positive and negative point vortices. Chavanis consid- 
ered a gas of point vortices with the same circulation. In 
this paper, we are interested in Chavanis's results. 

The system considered by Chavanis is closely related 
to the system of non-neutral plasmas in the Malmbcrg 
trap El El, HE M, SI SE Q ■ The dynamics of the 
non-neutral plasma in the Malmberg trap is described 
by the 2D Euler equation in a circular domain, using 
the guiding center approximation. Thus, if we are in- 



terested in non-neutral plasma in the experimental sit- 
uation, as an idealized model, we should consider a 
point vortex system in a circular domain. From experi- 
ments on non-neutral plasmas, many interesting proper- 
ties of vortex dynamics of the 2D Euler equation are now 
known El El, SE SH SI SE S3 : (!) Diocotron insta- 
bility (i.e., in other words, Kclvin-Hclmholtz instability), 
(2) Violent relaxation, (3) Slow decay, (4) Vortex crys- 
tals, and (5) Merger of vortices. As a theoretical aspect, 
recently the slow decay was numerically analyzed using 
the point vortex system [30l. l3ll . [Hj]. Although there is 
a difference in the boundary condition, these properties 
introduced here are common in the point vortex systems 
in an infinite plane. 

Back to the point vortex systems in an infinite plane, 
Chavanis derived serveral kinetic equations for the point 
vortex system in an infinite plane, in which the point vor- 
tices have the same circulation, and estimated interest- 
ing physical quantities, like the diffusion coefficient and 
the drift term, by using his kinetic equations EE EE E3 • 
In this paper, we develop a kinetic theory through the 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hier- 
archy [IE H, SE SI (for the text, See (HI). We compare 
our theory with that by Chavanis. The main purposes of 
this paper is as follows: (1) We derive a kinetic equation 
(i.e., the Landau equation), which is identical to that by 
Chavanis (Tol. E3- (2) Futhermore, we derive a kinetic 
equation (i.e., the Balescu-Lenard equation), which in- 
cludes more correlation, i.e., the collective effects. This 
is a new kinetic equation. (3) We show that for large N, 
this kinetic equation is reduced to the Landau equation. 

The organization of this paper is as follows. In § [TT1 
the equations of motion for a point vortex system in an 
infinite plane are presented. In § IIII1 the treatment of 
the BBGKY hierarchy is shown. Two key equations for 
time-evolution of the distribution function and the cor- 
relation function are derived. Using these equations, we 
derive the Vlasov equation, the Landau equation and the 
Balescu-Lenard equation for point vortex system in an 
infinite plane. It is shown that for large N, the Balescu- 
Lenard equation is reduced to the Landau equation. In 
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§ IIV1 we summarize the results of this paper and give 
future problems. 



II. EQUATIONS OF MOTION 

Let us consider a point vortex system, which consists 
of N point vortices with the same circulation 7 in an 
infinite plane. The hamiltonian of this system is given 
by 



2 N 



(1) 



where = (xi,yi). The equations of motion are written 
by using the hamiltonian: 



dxi dH dyi dH 
^ dt dyt ' ^ dt dxi 



(2) 



Thus, the velocities of the ith point vortex in the x- and 
y-directions are given by 



N 



E 



Xv) 



7 \ - (ri-rj)x 



(3) 
(4) 



It is convenient to rewrite the velocity in the following 
form. 



vt = E v ^' 



(5) 



where 



7 1 
2tt h-r^f 



J is the 2x2 symplectic matrix 

1 



J = 



-1 



(6) 



(7) 



z is the unit vector along the z-axis, and ry = — r^. 
This system have several conserved quantities: (1) 
Energy, i.e., H = E. (2) Angular impulse, i.e., 



I = lJ2iLi \ r i\ 2 - ( 3 ) Lm ear impulse, i.e., L = 7l]ili r * 
Now we define the vorticity w(r; t), 



N 



(8) 



and the stream function ip(r;t), 



^(r;i) = -^^ln|r-ri| 

Z7T * — ' 

1=1 



(9) 



Using the stream function, the velocity of the ith point 
vortex is given by 



Vj = -z x V^(r = Ti;t) 
= J- W(r = r<;t). 



(10) 



It is easily confirmed that the vorticity satisfies the two- 
dimensional Euler equation. 



~dt 



v • Vw = 0. 



(11) 



III. THE BBGKY HIERARCHY 

We define the distribution function of iV-point vortex 
systems: 



F = F(ri,r 2 , . . .,r N ;t). 



(12) 



The Liouvile equation for iV-point vortex systems in an 
infinite plane is given by 



where 



JY 



dt 



1 (dH d dH d 



(13) 



! = 1 



7 V dxi dyt dyi dx l 



(14) 



The Liouvillian C is rewritten as follows. 

i=l V yl 1 ' 



N 



(15) 



where Vj = and V y - = V, — Vj. Now we have used 
the fact that v(j — > i) = — v(i — > j). We define Cij as 

A? = - v (i *) ' V y 

= -v(j - i) • V; - v(t - j) ■ Vj. (16) 

Thus, the Liouvillian becomes 

JV 



i<j 



(17) 



We define the s-body reduced distribution function: 
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AP f 

/ s (ri,...,r s ) = — - — J dr s+1 dr s+2 . . .dr N F(r 1 , . . . ,r s ,r s+1 , . . . ,r N ). 



(18) 



r 



Carrying out a usual manipulation of the BBGKY hier- 
archy (See 29] in detail), we obtain the following time- 
evolution equations for the s-body reduced distribution 
function f s (r x , . . , , r s ): 



and 



dtfo = 0, 
d t f s (ri, . . .,r s ) 

£ijfs( v i, ■ ■ ■ ! r s) 



(19) 



i<j 



J dr s+1 £ iiS+ i/ s+ i(n, . . . ,r s+ i). (20) 



This is the BBGKY hierarchy for the s-body reduced dis- 
tribution functions. The time-evolution of the s-body re- 
duced distribution function is determined by the s-body 
reduced distribution function and the (s + l)-body re- 
duced distribution function. 
For s = 1, we have 



d t fi(ri) = J dr 2 £i 2 / 2 (ri, r 2 ) 



(21) 



Now we use the 2-body and 3-body correlation functions: 
52(ri,r 2 ) and g 3 (ri, r 2 , r 3 ). The 2- and 3-body correla- 
tion functions are related to the 2- and 3-body reduced 
distribution functions as follows: 

/ 2 (ri,r 2 ) = / 1 (r 1 )/ 1 (r 2 )+ 32 (r 1 ,r 2 ), (22) 
/ 3 (ri,r 2 ,r 3 ) = /i(ri)/i(r 2 )/i(r 3 ) + A(ri) 52 (r 2 , r 3 ) 
+/i(r 2 )g 2 (ri,r 3 ) + /i(r 3 )g 2 (ri,r 2 ) 
+ff3(ri,r 2 ,r 3 ). (23) 

Thus, the correlation function describes the deviation 
from the product of the reduced distribution functions. 
Inserting eq. (f2"2"|) into eq. (|21[) . we obtain 

<9 t /i(ri) = J dr 2 (£i2/i(ri)/i(r 2 ) + £i 2 # 2 (ri,r 2 )) . 

(24) 

Inserting eq. (f2"2"]l and eq. (|2"3")) into eq. (|2"0)) . the time- 
evolution equations of / 2 (ri,r 2 ) is obtained: 

9t/ 2 (ri,r 2 ) 
= £i2[/i(ri)/i(r 2 ) + g 2 (r 1 ,r 2 )} 

+ J dr 3 {Aal/iCri)/!^)/!^) 

+/l(ri)32(r 2 ,r 3 ) + /i(r 2 )g 2 (ri, r 3 ) 
+A(r 3 ) 52 (n,r 2 ) + 53 (n,r 2 ,r 3 )] + (1 & 2)}(25) 



Similarly, using eq. inserting eq. (f2"2"|) into eq. 1(23]). 
and manipulating the resulting equations, the time- 
evolution equation of g 2 (ri,r 2 ) is obtained: 

<9t5 2 (ri,r 2 ) 
= £i2/i(ri)/i(r 2 ) +£i2#2(ri,r 2 ) 

+ y^ r 3 {-Ci3/i(ri).92(r2,r3) + £23./i(r2).92(ri,r3) 
+(£i 3 + £ 2 3)[/i(i"3)5 2 (i'i,r 2 ) +5 3 (ri,r 2 ,r 3 )]}.(26) 

Equation (f2"oT) is slightly different from the recent result 
by Chavanis[17j|. who also developed a BBGKY hierarchy 
for the point vortex gas. The reason for these differences 
is unknown. 

Following [ItJ, we shall close the hiearchy of BBGKY 
equations by considering an expansion in powers of 1 /N 
for TV — > +oo. In the large AMimit, i.e., the hydrody- 
namic limit, we should preserve the total circulation T. 
Thus, the circulation 7 is 7 = T/N, where T — const. 
The order estimate of each function is / ~ 1, g 2 ~ 
l/N, 7 - l/N, and C %j ~ 1/AT. Thus, in eq. flSBJ), the 
first term is ~ l/N. The second term is ~ l/N 2 . The 
integral part is ~ l/N, since the integration over r 3 gives 
a A''- factor. The term including the term g 3 is ~ l/N 2 . 
In the following treatment, the function y 3 is omitted, 
since we cut the correlation, i.e., truncate a chain of the 
BBGKY hierarchy. 



A. Vlasov equation 

A mean field approximation is performed by neglecting 
the term of the correlation function in eq. (j2"4"| : 

dtfin) = J dr 2 £ 12 /(n)/(r 2 ). (27) 



Using eq. p^|) . we obtain 



at 



+ (vi) • Vi/(ri) = 0, 



(28) 



where 



(vi) 



dr 2 f(r 2 ){v x (2 l)i + v y (2 -► l)j). (29) 



This is a mean field equation for the point vortex sys- 
tem in an infinite plane. It is analogous to the Vlasov 
equation in plasma physics and in stellar dynamics. We 
should note that this equation, eq. I|28p. is nothing but 
the 2D Euler equation. 
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B. Landau equation 



(-(vi) ■ Vi - (v 2 ) ■ V 2 )s2(ri,r 2 ). 



(31) 



The next higher order approximation is started with Therefore, the correlation function is advected by (vi) 

and (V2). Equation (|3"Tjl is formally solved as 



eq. (I24p preserving the correlation function 

d t f(r 1 ) = J dr 2 (£i2/(r 1 )/(r 2 )+£ 1252 (r 1 ,r 2 )). (30) 

The right hand side of eq. (|30[) is upto the order 1/N 2 . 
For the correlation function gi 2 (ri,r 2 ), we approximate 
eq. ([21]) upto the order 1/N. 

9t52(ri, r 2 ) 
= £ 12 /(n)/(r 2 ) 



dr 3 (Ci 3 +£ 2 3)/(r 3 )32(ri,r 2 ) 
£i 2 /(ri)/(r 2 ) + 



Sf2(ri,r 2 ;t) 
C/i 2 (i)32(ri,r 2 ;0) 



/' 

Jo 



dr U 12 (T)C 12 f(r i; t- T)f(r 2 ;t - r), (32) 



where 
U 12 (t) 



exp 



dt'(vi)- Vr - / dt'(v 2 )-V03) 



Inserting eq. (|3"!2)) into eq. (|30p . we obtain 



fit/(ri) = J dr 2 C 12 f{r 1 )f(v 2 ) + J dr 2 C 12 U 12 {t)g 2 {v 1 ,Y 2] 0) 

+ J dv 2 dr C 12 U 12 (T)£. 12 f(r i; t - r)f(r 2 ;t - r). (34) 
In the right hand side of eq. (fM)) . the second term vanishes for large t (i.e., the correlation decays.), thus we have 

d t f(rx) = J dv 2 £i 2 /(ri)/(r 2 ) + ^ dr 2 J dr A2^2(r)r 12 /(ri;t - r)/(r 2 ;t- r). (35) 
We have to evaluate the following term: 

KconU/l = / dr2 / ^£ 12 [/ 12 (r)A2/(ri;t-T)/(r 2 ;i-T). (36) 

Using the approximation, 

/(r i; t - r)f(r 2 ;t - r) « (7 12 (-r)/(n; i)/(r 2 ; t), (37) 

we obtain 

^coiiU/} « / *2 J dr£i 2 U 12 (T)Ci2U 12 (-r)f(rr,t)f(vr,t) 

= Vi jdr 2 J dr v^v^t - T)V 12 f(r i; t)f(r 2 ;t). (38) 

The kinetic equation obtained here is 

+ (vx) • Vi/(r x ) = Vi • y^r, j dTV!(t)vi(t-T) • V 12 /(r i; i)/(r 2 ;i), (39) 

where Vi(t) is advected as vi(t — r) = £/i 2 (r)vi(i)/7i 2 (— r) and r^(t — r) = — /J" di' (Vj)(rj(t — — t'). This 
equation is analogous to the Landau equation in plasma physics and in stellar dynamics. This equation coincides with 
the result of [HI, 17 1. As shown in [l5j], this equation conserves the angular impulse and the linear impulse. If we use 
the Markovianization, i.e., extending the time integral to infinity, we obtain 



KcoiiU/} » Vi'fdea J drvi(i)vi(t-r).Vi 3 /(ri;t)/(r2;t). 



(40) 



However, it is not known whether the Markovianization gives long-time tail, i.e., the strong correlation. In par- 
is assured or not, since point vortex dynamics sometimes 
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ticular, in 31], it is shown that the diffusion process for 
the point vortex exhibits Levy flight. 

Chavanis estimated the relaxation time_£r e iax by us- 
ing the estimate of the diffusion coefficient [la] as i re iax ~ 
iV/ (In N)tn, where the dynamical time is to ~ (w) _1 ~ 
i? 2 /r, which is the time determined by the mean rota- 
tion time, and R is the size of the vortex. His estimate 
of [15J would be incorrect. In the kinetic theory, the N- 
dependence of t re i ax is determined by the TV-dependence 



of the collision term, i.e., IC 



(L) 
coll 



^reiax ~ Ntp. Recently Chavanis and Lemou used this 
estimate fill [l7[. This estimate is consistent with the nu- 



merical result by Kawahara and Nakanishi for the system 
in a circular domain 31 1 . 



C. Balescu-Lenard equation 

In this subsection, we derive a kinetic equation for 
point vortex systems in an infinite plane, which is analo- 



gous to the Balescu-Lenard equation in plasma physics. 
The starting point is the time-evolution equations of the 
one-body reduced distribution function f{v{) and the 
two-body correlation function g 2 {Yi,Y 2 ). 



dtf(n)= fdr 2 (£ 12 /(r 1 )/(r 2 )+A2ff 2 (ri ! r 2 )), (41) 



0(1/N). This gives and 



<9iff 2 (ri,r 2 ) 
= £i 2 /(ri)/(r 2 ) + A 2 9 2 (ri,r 2 ) 

+ J dr 3 {£13/(1*1)32 (r 2 ,r 3 ) + £ 23 /(r 2 )# 2 (n, r 3 ) 

+(A3 + £ 23 )[/(r 3 )s 2 (ri,r 2 )]}. (42) 

The two-body correlation function is formally solved as 



52 (ri,r 2 ;t) = / dr U 12 (t) {C 12 U 12 {-r)f{T 1 )f{v 2 ) + C 12 U 12 (-r)g 2 (r u r 2 ) 
Jo 

+ J dr 3 (£i 3 C/i 2 (-r)/(ri) 52 (r 2 ,r 3 )+£ 23 C/i 2 (-r)/(r 2 ) ff2 (ri !r3 ))} (43) 

I 



The kinetic equation is formally obtained as 

^ + (vi)-Vi/i= / dr 2 C l2 g 2 {v u v 2 -t). (44) 

If, as done for plasma systems in [2^], we set 
5 2 (ri,r 2 ;t) = .g 2 (ri - r 2 ;t) 

= Jdkexp[ik- (n -r 2 )]g 2 (k;t), (45) 



the right hand side of eq. ([44|) . i.e., the collision term, 
vanishes. 

J dr 2 £i25 2 (i"i,i" 2 ;0 
= - Jdr 2 (V 1 F(r 1 -r 2 )) T - J-V l52 (ri,r 2 ) 

J dr 2 J dk exp[ik • (r x - r 2 )]V(k)(ik) T ■ JVi 
J dk' exp[ik / -(n-r 2 )]5 2 (k';t) 



(2tt) S 



dkV(k) k T ■ J -kg 2 {-k;t) 



= 0, 



(46) 



since k T ■ J ■ k = 0. Here "A T " means the transpose of 
the vector A. In eq. (|4"B")) . the homogeneity is assumed. 



The above result shows that inhomogeneity is important 
for point vortex systems. To not make the collision term 
vanish, we change the definition of the Fourier transform 
of <7 2 (ri,r 2 ;t). Therefore, we set 

52(1-1, r 2 ;t) 

= J dk! J dk 2 exp[iki • n + ik 2 • r 2 ]g 2 (ki, k 2 ; /()47) 

Ci 2 can be rewritten in the form 

C 12 = -(V 1 V(r 1 -r 2 )) T ■ J- V12, (48) 
7 



where 



V(rx - r 2 ) = - — In |r x - r 2 | 

Z7T 



(49) 



Now we consider the Fourier transform of the function 
V(r), where r = ri — r 2 : 



V(r) = I dk V(k)e l - 
1 



V(k) 



drV(r)e~ 



(50) 



The Fourier transform of V^(r) is evaluated as follows: 
V(k) 



6 



7 



(2tt) 3 

,/2 



7 



(2tt) 3 



dr In |r| e - lk r 

2tt 



limiting value of this function is zero. Alternatively, we 
insert a convergence factor: 



rdr d6 lnre^ fcrcos(f,) 

n Jo 



. , , dr rlnrJo(fcr) 
(2tt) 2 



(2^Ufc lnrJl(fcr) 



We have to evaluate the following limit: 



lim — lnrJi(fcr) 

r — >oo /c 



lim 

e^+0 



dr rlrtrJo(kr)e 



(53) 



dr Ji(fcr) bl) As a result ' we have 



V(k) = 



(52) 



7 2 1 
(2tt) 2 fc 2 ' 



(54) 



In fact, for large r, this function oscillates with ampli- The dependence of V(k) ~ 1/fc 2 is a typical behavior of 
fying its absolute value. Therefore, we suppose that the Coulomb systems. 



J 



Now we evaluate each term in the right hand side of eq. ((43 



(A) 



dTU 12 {T)C 12 U 12 {-T)flh 

f drvi(t-r)- V12/1/2 
Jo 



1 Jo 



dr / dkexp 



ik- (n - r 2 ) -ik- / dt'{vi)+ik- dt'(v 2 ) 



V(k)(ik) T ■ J-Vi 2 /i/ 2 . (55) 



/ dr C/i 2 (T)£i 2 i7i2(-T)52(ri,r2;t) 
Jo 

-/ dr vi(i) • Vi2S , 2(ri,r 2 ;t) 
Jo 



— dr dkexp 
1 Jo 



ik-(r 1 -r 2 )-ik- dt'(vi) + ik • / dt'{v 2 ) 



V(k)(ik) T ■ J ■ V± 2 g 2 (ri, r 2 ; t). (56) 



(C) 



J dTU 12 (r) J dr 3 £ 13 U 123 {~T)f 1 g 2 (r 2 ,r 3 ;t) 



dR / dR' / dr dk dk' 



(2tt) 2 7 



x exp 



ik • (n - R') + zk' ■ (r 2 - R) - ik • J ((vi> - (V'))dt' 



xf(A;)((ik) T • J-V 1 / 1 ) ff2 (R,R';t). 



(57) 



(-D) = J dTU l2 {r) J dv 3 C 23 U l23 {-T)f 2 g 2 {v u v 3 ;t) 



(2tt) 2 7 



dR / dR' / dr dk dk' 



x exp 



ik' • (n - R) + ik ■ (r 2 - R') — ik - J «v 2 ) 



<V»dt' 



xT/(fc)((»k) T • J- V 2 / 2 ) ff2 (R,R';t)- 



(58) 



7 



The order estimate of these terms is as follows. 



(A) ~ — , IB) ~ -At, (C) ~ — , (D) ~ — . 



Therefore, we can neglect the term of (B). 

Then we obtain the following integral equation: 



where 



and 



<7 2 (ri,r 2 ;t) = g(n,r 2 ;i)+ J dR J dR' K({r 1: r 2 }, {R, R'})g 2 (R, R'; t), 



1 

Uo 



J dr J dkexp ik ■ (ri - r 2 ) - ik • J dt'(vx) + ik ■ J dt'{v 2 ) 



xV(k)(ik) T - J-VuA/j, 



K({r 1; r 2 },{R,R'}) 

dr J dkj dk' exp ik • (n - R) + ik' • (r 2 - R) - ik • J ((vi) - (V'))<2*' 



(27T) 

xV(k)({k) T ■ J ■ Vi/; 



so / dr [ dk [ dk' exp 

7r) 2 7 7o 7 7 



(2 

xf(fc)(ik) T - J-V 2 / 2 . 



ik' • (ri - R) + ik • (r 2 - R') 



ik 



((v 2 ) - <V'»df 



r 



(59) 



(60) 



(61) 



(62) 



The function K({ri, r 2 }, {R, R'}) is called the integral 
kernel of the integral equation. This integral equation 
takes the form of the Fredholm integral equation of the 
second kind. Thus, how to solve it is known -33]. Now 
for brevity, we set x = (ri,r 2 ) and y = (R, R'). The 
integral equation, which should be solved, is 

0»(x;t) = q(*;t) + \ J dyX(x,y). 92 (y;t). (63) 

If the required conditions are satisfied, this integral equa- 
tion is solved as 



The function H(x, y; A) is called the resolvent. Ao is cho- 
sen to make the series convergent. For large N, the ker- 
nel if(x, y) is ~ 1/N. Thus, for large N, we can take as 
Ao = 1. Therefore, if the kernel is bounded, for large N, 
we get a convergent series. The resolvent is given by 



S(x,y; A) = 



^(x,y;A) 
D(X) ■ 



(65) 



02(x;t) = q(x;t)+ / dy H(x,y; X )q(y;t). (64) where 



D(X) = 1 — X ds K(s,s) 



(-X)p 



and 



2! 



0(x,y;A) = 



ds< 



tf(x,y)-A JdsK^* 
, (-A) p 



dsids 2 K 



■ ds p K 



dsi ■ ■ ■ ds p K 



si s 2 
si s 2 



x si 

y s 2 



(66) 



% [ [d Sl ds 2 K(* 81 82 

21 J J \ y s 2 s 2 



X S! 

y s 2 



(67) 
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and 



K 



Si s 2 • ■ • s p 

ti t2 ■ ■ ' t p 



K(s u ti) K(s u t 2 ) ■■■ K( Sl ,t p ) 
K{s 2 ,t 1 ) K{s 2 ,t 2 ) ••• K(s 2 ,t p ) 

#(s p ,ti) K(s p ,t 2 ) ••• K(s 



(68) 



Thus, a derived kinetic equation is 



dt 



i f 

- + (vi)-Vi/i= dr 2 C 12 g 2 (n,r 2 ;t) 7 



where 



<fc(ri,r 2 ;t) = <j(ri,r 2 ;t) + / dR / dR' S({ri, r 2 }, {R R'}; A ) <?(R, R'; i). 



(69) 



(70) 



This kinetic equation for point vortex system is an ana- 
logue of the Balescu-Lenard equation in plasma physics. 
Unlike the Balescu-Lenard equation for plasmas, the col- 
lision term obtained does not have the dielectric func- 
tion, but have additional terms compared with the Lan- 
dau collision term. We note an important point, which 
manifests a difference between point vortex systems and 
plasma systems. The kernel has the factor 1/N, which 
comes from the circulation 7 = T/N. Thus, the resolvent 
is expanded in terms of the 1/iV-factor. The integral is 
just the Fourier transform. Therefore, the integral is or- 
der of 1. In eq. the second term of the right hand 
side is smaller than the first term as the order 0(1/N). 
Therefore, for large N, we obtain 

g 2 {ri,r 2 ;t) q(n,r 2 ;t). (71) 

With this result, the collision term is reduced to the 
Landau collision term, i.e., eq. (|38|) . Therefore, for point 
vortex systems, the Balescu-Lenard collision term is 
reduced to the Landau collision term, because of large 
N-eSect, i.e., the absence of the dielectric function. This 
is a big difference between point vortex systems and 
plasma systems. The screening effect does not appear 
in point vortex systems, while the screening effect is 
important for the Balescu-Lenard equation for plasmas. 
On the Balescu-Lenard collision term in plasma physics, 
see [1^] for comparison. 

In the above paragraph, we showed how to solve the 
integral equation eq. (j6T?l) formally. But the required con- 
ditions have not yet been checked. Let us look at the 
required conditions. First, to show the operator T 

Tg 2 = J dy *r(x,y) 52 (y), (72) 

is bounded, is important. The integral equation is sym- 
bolically given by 

(1 - XT)g 2 = q. (73) 



The formal solution is obtained by the expansion. 

.g 2 = ^(AT)" g . (74) 

The Fredholm theory for integral equation basically uses 
the boundedness of the operator and the expansion, i.e., 
eq. (|74|) . This expansion is assured by the boundedness of 
the operator T. The expansions of eqs. (f6"6")) and (|6"T|) are 
also due to the boundedness of the operator T. As the 
second point, unlike usual Fredholm integral equations, 
the integral domain in eq. (|6li|) is infinite. If the integral 
domain is finite and the kernel is bounded, the bounded- 
ness of operators is easily shown. For our case, since the 
integral domain is infinite, we should treat the operator 
T carefully These two points, i.e., the boundedness of 
the operator T and the infinite integral domain, should 
be checked and be treated in a rigorous way. However, 
in this paper, we do not pursue a rigorous discussion. 
These problems are reserved for mathematical physicists. 



IV. CONCLUDING REMARK 

We have derived a kinetic equation for point vortex sys- 
tems in an infinite plane. The kinetic equations derived 
are analogues of the Landau equation and the Balescu- 
Lenard equation. Equation (|39[) coincides with the re- 
sult of Chavanispjl [lTj]. Equation (f39|) possesses sev- 
eral interesting properties, which were shown in [l5l |. 
The Balescu-Lenard equation (|69[) is a new kinetic equa- 
tion. Furthermore, we have shown that for large N, the 
Balescu-Lenard equation, i.e., eq. (|69|) is reduced to the 
Landau equation, i.e., eq. (|39p . Therefore, we can con- 
clude that for point vortex systems in an infinite plane, 
without symmetrical restriction (such as axisymmetric 
and unidirectional flows), the most generalized kinetic 
equation is eq. (f3"9")l . 

The following point would be the interesting point of 
the derived kinetic equation, i.e., eq. (|39p : The interac- 
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tion among point vortices is long range, i.e., logarithmic. 
In addition, the derived kinetic equation is analogous to 
the Landau equation for 3D plasmas. However, the col- 
lision term for point vortex systems may not diverge. In 
[l7l ]. for the axisymmetric case, the collision term does 
not diverge. This is a symptom of the non-divergence 
of the collision term for the Landau equation. The rea- 
son of this is that the difference between the integration 
of the collision term, i.e., eq. (|38|) . and that of the Lan- 
dau collision term for 3D plasmas. The former has the 
integration with respect to r 2 (i.e., position), while the 
latter has the integration with respect to r 2 and v 2 (i.e., 
velocity) . 

For the point vortex systems, the expression of the 
energ y sp ectrum was derived for the system in an infinite 
plane[34j and for the system in a circular domain (35j. 
The ener gy s pectrum is closely related to the diffusion 
coefficient [361] . We will be able to compare the theory of 
[3(| with the kinetic theory in this paper. 

Another interesting point is the following: our kinetic 
theory is not directly connected to Onsager's tempera- 
ture. As shown in [371 ] for two-sign point vortex systems, 
Onsager's temperature affects nonequilibrium properties, 
i.e., decaying process. To find this connection leads to 
understand nonequilibirum properties, i.e., classification 
of nonequilibrium processes. 

The most interesting problem, which we would like to 
attack with eq. ([39]) . is a decaying property of vortex 
crystals in non-neutral plasmas. Unfortunately, eq. (|39|) 
is for the system in an infinite plane, not for the system 
in a circular domain. Thus, it is not for an experimen- 
tal situation. But, eq. (|39[) surely captures the nature of 



phenomena for the system in a circular domain in some 
extent. Many experimental results show that the vortex 
crystals are quasi-stationary states. To analyze quasi- 
stationary states, recent advances for long-range inter- 
action systems would be some hints for us, such as a 
study of the Hamiltonian mean field (HMF) model. The 
HMF exhibits a slow decay, in which the state is stuck 
in a quasi-stationary state, as well as in point vortex sys- 
tems. For the HMF model, the slow decay is analyzed by 
the Vlasov equation. The quasi-stationary state is very 
near to a stable solution of the Vlasov equation. The 
estimate of some quantities, i.e., the algebraic decay and 
the tail of the velocity distribution function etc., is car- 
ried out [H|], and is tested by a numerical simulation [39j. 
Their analysis would be useful for our problems. How- 
ever, to this end, we should know the bahavior of the 2D 
Euler equation, instead that of the Vlasov equation for 
ususal kinetic theory of particle systems. 
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